TTP97-05, 



hep-ph/ 9703291 



Two-loop QCD corrections to semileptonic b decays at maximal 

recoil 



OS 



43 

Oh! 



Andrzej Czarnecki and Kirill Melnikov 

Institut fur Theoretische Teilchenphysik, Universitdt Karlsruhe, 
D-76128 Karlsruhe, Germany 
(March 1997) 



5h 

oj ; Abstract 



We present a complete 0{a 2 ) correction to the differential width of the in- 
clusive semileptonic decay b — > c/i/; at the kinematical point of vanishing 
invariant mass of the leptons, q 2 = 0. Together with the recently computed 
correction at the upper boundary of the lepton invariant mass spec- 
OS ■ trum, this new information permits an estimate of the 0{at 2 s ) effect in the 

■ total inclusive semileptonic decay width b — ► clv[. We argue that the non- 

BLM part of the 0(a 2 ) correction gives at most 1% correction to the inclusive 
semileptonic decay width b — > clv\. This significantly improves the credibility 
of extracting |V^,| from the inclusive semileptonic decays of the 6-hadrons. 



Semileptonic decays of the b quarks provide the best opportunity to determine \V c b\, 
a parameter of the Cabibbo-Kobayashi-Maskawa (CKM) matrix and a fundamental input 



parameter of the standard model. The current experimental limit |IJ] 

\V cb \ = 0.036 to 0.046 (90% CL) (1) 

H , 

is based on measurements of the beauty hadron decays produced at the T(4S I ) resonance 
(by ARGUS and by CLEO II) and in Z-boson decays (by the four experiments at LEP). 
In the future large samples of the 6-hadrons collected at 5-factories (at SLAC and KEK) 
and at the hadron colliders will increase the statistical accuracy to a few percent level. To 
fully exploit the anticipated experimental improvement, the theoretical description of the b 
decay must be known with comparable precision. 

There are two methods of extracting the value of \V c b\, based on measurements of the 
exclusive decay B — > D*lui and of the inclusive semileptonic decay width of 6-hadrons 
T s \. These two methods rely on very different theoretical considerations and experimen- 
tal procedures and complement each other. Their merits and theoretical uncertainties are 
summarized e.g. in Ref. [@-|]]- One of the major sources of the theoretical error are the per- 
turbative QCD corrections at the two loop level. For the exclusive decays at the zero recoil 
point these corrections have recently been calculated p. This has significantly improved 
the accuracy of the theoretical prediction for the exclusive method. 
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In the case of the inclusive semileptonic decay width of the 6-hadrons T s i, the only known 
effects beyond one loop are those associated with the running of the strong coupling constant 
|||7]]. They are obtained by computing massless quark effects (fig. 1) and then replacing the 
number of light flavors Nl by the combination in which it enters the one-loop /3-function 
Nl — 33/2. These so-called BLM corrections || are expected to dominate the two- loop 
result, however only a full calculation of the remaining diagrams will put this statement on 
a firm foundation. 

Technically, the correction to the semileptonic decay width T 8 \ is obtained by fixing 
the invariant mass of the leptons q 2 and computing the differential width dr s i/dg 2 with 
desired accuracy. Integrating over q 2 within kinematical boundaries, one gets the inclusive 
semileptonic decay width of b — > duf. 

{m b -m c ) 2 

r,= / tf^. (2) 

^ 

Going beyond the BLM approximation and computing complete 0(a 2 ) corrections re- 
mains a daunting task at present. In comparison with the zero recoil calculation the main 
difficulties are: an additional kinematical variable describing the invariant mass of the lep- 
tons (q 2 ) and the presence of the real radiation of one and two gluons. 

To circumvent these difficulties, we propose to estimate the deviations from the BLM 
predictions by performing complete 0(a 2 ) calcualtions for dr s i/dg 2 at two boundaries of 
integration in eq. (fj). 

In fact, one of these calculations has already been done in ref. H where 0(a 2 ) corrections 
to the transition b — > clvi were calculated at the zero recoil limit. Since in this limit the 
radiation of real gluons is absent, the results of || provide 0(a 2 ) correction to dr sl /dg 2 at 
?L = Oft ~m c ) 2 . 

The purpose of this paper is to present a calculation of the 0(a 2 ) corrections at g^ in = 
which is the other boundary for the invariant mass of the leptons. With both boundary 
points known we can estimate the deviation of the 0(a 2 ) corrections to the total inclusive 
semileptonic decay width of the 6-quark r s i from the BLM prediction. 

Taking q 2 = limit is important for the feasibility of this calculation. In this case the 
calculation of real radiation of one and two gluons is considerably simplified. 

The reason why the real radiation at order 0(a 2 ) is difficult to calculate is that the 
particle in the initial state (the decaying b quark) carries a color charge and therefore can 
radiate. It is the presence of the massive propagators of this particle which makes the 
integrations over the phase space very tough. For this reason even the QED corrections to 
such well studied processes as the muon decay remain unknown at the two-loop level. The 
kinematical configuration in which q 2 = and the quark in the final state is massive is the 
first case where the complete evaluation of the real radiation in the decay of a fermion turns 
out possible. Below we sketch the basic ideas of our approach; the technical details will be 
presented elsewhere. 

The idea which permitted us to calculate the contribution due to the real radiation of 
one and two gluons is (qualitatively speaking) the expansion in the velocity of the final 
quark. Indeed, in the limit m c — > nib the charm quark in the final state is a slowly moving 
particle, with spatial components of its momentum of the order of — m c , much smaller 
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than its mass. The four momenta of gluons and of leptons (for q 2 = 0) are also of the 
order of rnj, — m c . It turns out that by a proper choice of the phase space variables one can 
systematically expand the amplitudes and the phase space in terms of 5 = [m^—m^jm^ <C 1. 

Some examples of the diagrams which contribute to the QCD corrections to the semilep- 
tonic decay of the b quark are shown in fig. 2. Not shown are several other virtual corrections 
as well as diagrams obtained by permuting the gluon couplings to the quark line or by cross- 
ing the external gluon lines. In total there are about 80 Feynman diagrams which have to 
be evaluated. 

We do not include the diagrams with three c-quarks in the final state in our analysis. 
Since 3m c is only marginally smaller than m&, the contribution of such diagrams is strongly 
suppressed. 

We parametrize the expansion using the variable 6=1 — m c /mb- In the first two 
nonvanishing orders (S 3 and 5 A ) only virtual corrections contribute (e.g. fig. 2a,b). The 
following two terms receive in addition contributions from diagrams with one loop and one 
real gluon emission (like in fig. 2c, d), as well as from diagrams with two gluons resulting 
from a decay of a virtual gluon (fig. 2f). Only in the order S 7 the contributions of a double 
gluon emission from the quark line show up (fig. 2e). This hierarchy can be traced back to 
the fact (evident in physical gauges) that the interaction of the slowly moving quarks with 
real gluons is proportional to the three velocity of the former. 

In case of two-loop virtual corrections as well as in the emission of two real gluons the 
expansion in 5 means a Taylor expansion in the small external momenta of the leptons and 
gluons. Such an expansion does not lead to any spurious ultraviolet or infrared divergences. 
The situation is different in the case of the single gluon radiation in diagrams where there 
is in addition one virtual loop (fig. 2c, d). There a naive Taylor expansion in the external 
gluon momentum leads to artificial infrared divergences which correspond to the on-shell 
logarithmic singularities of the one-loop diagrams. Therefore a more sophisticated approach 
is needed and the recently developed method of "eikonal expansions" P,|T0[ is used. 

To present our result we write the differential semileptonic decay width of the decay 
b — > civ at q 2 = as 



dr sl 



dq 2 



A B o m + -C F A 1+ (^) 2 CVA 2 

7T \ 7T / 



(3) 



where T = I Kfe| 2 and A Bor n,i,2 describe the m c /m b dependence in various orders in the 
strong coupling constant. 

.2 /^2n3 



Both Aboi-ii — (1 — m c/ m b) an d Ai are known in a closed analytical form [|IT| , |I"2"f . A 2 is 
the main result of the present paper. For the purpose of presentation we divide it up into 
four contributions according to the color factors: 

A 2 = 5 3 i(C F -C A /2)A F + C A A A + T R N L A L + T R A H ] (4) 

The last term, Ah, describes the contributions of the massive b and c quark loops. Top 
quark contribution is suppressed by a factor ~ ml/m 2 and has been neglected. 

For the SU(3) group the color factors are C A = 3, Cp = 4/3, Tr = 1/2. Nl = 3 is the 
number of the quark flavors whose masses have been neglected (u, d, and s). 
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We have computed the expansion coefficients of Af,a,l,h up to 5 8 , which for the phys- 
ical value of the charm and bottom masses gives an estimated accuracy of our numerical 
predictions better than 1% (for 5 = 1 — m c /m b rs 0.7). 

In the present paper we list the analytical results only up to 5 4 , while the numerical 
evaluation is done using the expansions up to 5 8 . Using the pole mass of the b and c quarks 
and expressing the one-loop corrections in terms of oijjgijnl) we find 

355 2 2 /89 2 

A A = h -7T + 61 7T 2 

36 3 V 8 

_ 2 / 2422517 1708, . 44, 2 . 8 257 2 \ 

+6 2 1 \n(25) n 2 (2<5) + -ci H 7T 

V 32400 45 K J 9 v ; 9 90 / 

r3 /2956607 854, . 22, 2 . 4 307 2 

+5 (^4800- - 45- ln(25) + T ln (25) " 9 Cl - 180" 
a( 5789957 4663, 2, 2/an 4 412 2 



+5 4 + ln(25) + - In 2 (26) + — c x + 

V 1323000 4725 V ; 5 v ; 45 1575 

23 8 8 2 c/23 , 2 \ e2 /1697 8, 22 359 



- e ■ f + r + Ht - 4C2 - 4 " J + * - 3 ln(M) + "T 2 + 135 

V 720 3 v ; 15 270 

r4 /4957991 1460, , n _, 16, 2 . orA 2 139 2 
+<T ln(25) H In 2 (25) + -c 2 vr 2 

V 396900 189 v ; 9 v ; 7 600 

14 . ,2/82217 544, , o _, 16, 2/nr , 16 2 

A L = 5 + 5 2 [ ln(2<5) + — In 2 (25) vr 2 

9 V 4050 45 v ; 9 v ; 27 

r o / 103667 272, , ne . 8, 2/an 8 2 \ 

+5 3 + \n(25) - - \n 2 (25) + — vr 2 

V 8100 45 K J 9 V ; 27 / 

r4 / 1322183 2404, . 8 , 2/or . 8 2 

+ * (l9^-l^ ln(25) + 45 ln (25) -135- 

A fl = - ^ + 5 (-74 + 8. 2 ) + 5 2 (^- ^ 
9 3 V / V 81 27 

. 3 / 33883 32, . 136 2 \ , 4 /3754 154 2 
+5 6 ln(25) H vr 2 + <T vr 2 

V 810 9 V ; 27 J V 405 135 



-7T 



(5) 



with ci = f Cs - vr 2 ln(25) and c 2 = f Cs - vr 2 In 2. 

We now turn to the numerical analysis of our result. Here the issue of numerical values 
for the quark masses becomes important. It is safe to assume that the pole mass of the 
6-quark lies between 4.6 and 5.1 GeV. The mass of the c quark is determined by m& — m c , 
obtained from the HQET calculations 0,||,0[J. We use m b - m c « 3.45 ± 0.10 GeV where 
the error bar is rather conservative. 

Accordingly, the numerical value of 5 changes within the range of 0.65 — 0.77. The 
numerical values for the function A 2 become: 

A 2 = -6.03, - 7.45(4), - 8.96, (6) 

for 5 = 0.65, 0.7, 0.75 respectively. 

The error estimate, shown for the central value of 5 = 0.7, is obtained by multiplying the 
last computed term by 3, which corresponds roughly to summing up the remainder of the 
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series in 5 assuming constant coefficients. This procedure overestimates the error because 
the coefficients in fact decrease (there is at most a logarithmic divergence at 5 = 1 caused 
by neglected diagrams with three real c quarks in the final state). 

Taken literally, the 0(a 2 ) corrections are quite large. However, as we will show below, 
the bulk of them is due to the BLM corrections. 

The BLM prediction with 4 light flavors of quarks gives the following results: 

Af LM = -5*A L T R (§-*) = "6-54, - 8.15(6), - 9.87, (7) 

for 5 = 0.65, 0.7, 0.75 respectively. 

By comparing the numbers in eq. (0) with those in eq. @ we conclude that the 
BLM correction accounts for most of the effect. We estimate the residual correction 
by subtracting the BLM piece from the exact correction. We get a residual correction 
(0.51, 0.7, 0.91)C_p(a; s /7r) 2 , which, using a s (m b ) = 0.23, gives numerically 0.5, 0.7, 0.8% 
correction relative to the Born rate for 5 = 0.65, 0.7, 0.75. 

Therefore, we arrive at the conclusion that at the lower boundary of the invariant masses 
of leptons g^ in = 0, the BLM piece of the 0{a 2 ) correction represents the complete result 
with an excellent accuracy. The remaining correction does not exceed the value of 1% even 
accounting for an uncertainty in input parameters. 

Finally, we would like to estimate the 0(a 2 ) radiative corrections to the total semileptonic 
decay width of the h quark. In the BLM approximation such corrections have been calculated 
in ref . . Therefore, we are only interested in the deviations from the BLM approximation. 

Our estimate of the non-BLM corrections to the inclusive width is based on the expec- 
tation that the largest deviation from BLM should occur at the maximal recoil limit, i.e. at 
q 2 = 0. To clarify this point, we note that the results of ref. || imply that at zero recoil 
limit ((/^ax = ( m b ~ m c) 2 ) the deviation of the exact result from the BLM approximation are 
very small. On the other hand, the results of this paper show that at q 2 = the non-BLM 
part of the correction grows with the decrease of the c-quark mass, i.e. with the increase in 
the phase space available for real gluon radiation. If one fixes the value of the c-quark mass, 
but varies instead the invariant mass of leptons q 2 , the strongest emission of real gluons will 
occur at the maximal recoil point of the spectrum, at q 2 = 0. It is for this reason that we 
expect the largest discrepancy between the BLM prediction and the full correction at the 
lower end of the q 2 distribution, for q 2 = 0. 

Turning to the estimate itself, from ref. [f| we know that at g^ ax = (m^ — m c ) 2 (zero recoil 
limit) the non-BLM correction to the differential width relative to the Born value is of the 
order of —0.1%. On the other end of the lepton invariant mass distribution the result of this 
paper implies a slightly larger, but also tiny deviation below 1%. We note that the change 
of sign of the non-BLM corrections cancels part of their impact on the total width. Taking 
the absolute value of the larger of the corrections at the boundaries as an upper bound we 
conclude that the non-BLM piece of the 0(a 2 ) corrections to the total semileptonic decay 
width b — > clvi should not exceed the value of 1%. 

The value of second order correction to the inclusive width depends on the adopted 
definition of the quark masses. Our result is presented in terms of the pole masses, which 
is a convenient choice for the corrections not associated with the running of the coupling 
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constant. It was argued in JT3] that such parametrization leads to small higher-order non- 
BLM corrections. Our result confirms this expectation. 

It is fair to say at this point that our estimate of the non-BLM piece of the corrections 
to the total inclusive semileptonic decay width based on the two boundary values can not 
be considered as a rigorous proof. Keeping in mind that the complete calculation of the two 
loop QCD corrections to the total decay width remains a very difficult task, a calculation 
of these corrections at some intermediate point qf nt for the differential inclusive semileptonic 
decay width of the 6-quark is highly desirable. If such a calculation confirms that the non 
BLM piece of the correction remains within the range set by its value on two boundaries, 
our estimate for the correction to the total semileptonic decay width of the b quark will be 
on a very safe ground. 

We are grateful to K. G. Chetyrkin, M. Neubert and N. G. Uraltsev for discussions and 
advice in the course of this calculation. We would like to thank Prof. J. H. Kiihn for his 
interest in this work and support. This research has been supported by BMBF 057KA92P 
and by Graduiertenkolleg "Teilchenphysik" at the University of Karlsruhe. 
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FIGURES 




(a) (b) 

FIG. 1. Diagrams involving a light quark loop (a) or real pair emission (b). Symbols <S> mark 
places where the virtual W boson can possibly couple to the quark line. 



(a) (b) 




(e) (f) 



FIG. 2. Examples of the two-loop gluonic QCD corrections to the decay b — > clv\ at zero recoil; 
a,b: virtual corrections; c,d: single gluon emission; e,f: emission of two gluons. Symbols ® mark 
places where the virtual W boson can possibly couple to the quark line. The left hand side diagrams 
are QED-like, while the right hand side ones are purely nonabelian. 
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